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Using a sample of 106 million $w(3686) decays, the branching fractions of w(3686) — 
yxXco, (3686) — yxe1, and (3686) — yxe2 are measured with improved precision to be 
(9.389 +0.014+0.332) %, (9.905 -£ 0.011 40.353) 96, and (9.621+0.013+0.272) %, respectively, where 
the first uncertainties are statistical and the second ones are systematic. The product branching 
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fractions of (3686) + Yyxa,xa > yJ/w and (3686) > vyx«e, xe» — *yJ/w and the branching 
fractions of x«i > yJ/w and xe» > *J/1» are also presented. 


PACS numbers: 13.20Gd, 13.40Hq, 14.40Pq 


I. INTRODUCTION 


The discovery of the J/«» in 1974 and soon thereafter 
of the charmonium family convinced physicists of the re- 
ality of the quark model np Since then, measurements 
of the masses and widths of the charmonium family and 
their hadronic and radiative transition branching frac- 
tions have become more precise. The spectrum of bound 
charmonium states is important for the understanding of 
Quantum Chromodynamics (QCD) in the perturbative 
and non-perturbative regions B. 

For charmonium states that are above the ground state 
but below threshold for strong decay into heavy flavored 
mesons, like the 15(3686), electromagnetic decays are im- 
portant decay modes. The first charmonium states dis- 
covered after the J/y and (3686) were the yc; (J = 0, 
1, and 2) states, which were found in radiative transitions 
of the (3686) [3,14]. These states, which are the triplet 
1P states of the cc system, had been theoretically pre- 
dicted E! [6] along with the suggestion that they could be 
produced by F1 transitions from the 1»(3686) resonance. 

Radiative transitions are sensitive to the inner struc- 
ture of hadrons, and experimental progress and theoret- 
ical progress are important for understanding this struc- 
ture. The development of theoretical models is also im- 
portant for predicting the properties of missing char- 
monium states, in order to help untangle charmonium 
states above the open-charm threshold from the myste- 
rious XYZ states i. Much information on radiative 
transitions of charmonium can be found in Ref. i, and 
a recent summary of theoretical predictions for radiative 
transitions of charmonium states and comparisons with 
experiment may be found in Ref. [s]. 

The branching fractions of 1/(3686) > yve7 were mea- 
sured most recently by CLEO in 2004 with a sample of 
1.6 M (3686) decays [9]. 'The Crystal Ball it CLEO 
values, and the Particle Data Group (PDG) [7] averages 
are given in Table [I] 


TABLE I. Crystal Ball and CLEO [9] (3686) — yer 
branching fractions and average values from the PDG [1]. 


Po Gn 
8 


Decay 

~(3686) > yxco| 9.92: 0.5 £0 9.22 + 0.11 + 0.46 | 9.2 + 0.4 
(3686) > yxcı | 9.0 + 0.5 +0.7 |9.07 +0.11 +0.54| 8.9 + 0.5 
~(3686) > yxc2| 8.0 + 0.5 +0.7 |9.33 +0.14 + 0.61 | 8.8 + 0.5 


BESIII has the world’s largest sample of 1»(3686) de- 
cays and has made precision measurements of many 
(3686) branching fractions, including (3686) — 
nv^7-J/V, along with J/y > I (L = e,u) [U], 
v(3686) — «?J/v and qJ/v [12], (3686) — 7h, (13, 


[14], and the product branching fractions B(2)(3686) > 
yXcJ) X B(xey > yI/W) [15] iid using exclusive XeJ — 
yJ/w decays. It is important that the (3686) > yxe7 
and (3686) — yne branching fractions be measured as 
well. Improved precision on these is necessary because 
they are often used in the determination of yey and me 
branching fractions via the product branching fractions. 
However, it is to be noted that systematic uncertainties 
dominate the measurements summarized in Table [I| so 
to improve on their results, it is necessary to reduce the 
systematic uncertainties. 

In this paper, we analyze (3686) inclusive radiative 
decays and report the measurement of the (3686) — 
YXeJ branching fractions. The product branching frac- 
tions B(wW(3686) — yver) x Blver — »yJ/v) are also 
measured, and the xez — yJ/w branching fractions 
are determined. This analysis is based on the w(3686) 
event sample taken in 2009 of 106 million events, deter- 
mined from the number of hadronic decays as described 
in Ref. [17], the corresponding continuum sample with 
integrated luminosity of 44 pb~! at V/s = 3.65 GeV [17], 
and a 106 million (3686) inclusive Monte Carlo (MC) 
sample. 

The paper is organized as follows: In Section II, the 
BESIII detector and inclusive (3686) MC simulation 
are described. In Section III, the selections of inclusive 
(3686) + yX events and 7°’s are described and com- 
parisons of inclusive (3686) data and MC sample dis- 
tributions are made. Section IV presents the inclusive 
photon energy distributions, while Section V details the 
selection of exclusive (3686) + ^x.; events. Sections 
VI and VII describe the fitting of the photon energy dis- 
tributions and the determination of the branching frac- 
tions, respectively. Section VIII presents the systematic 
uncertainties, and Sections IX and X give the results and 
summary, respectively. 


II. BESIII AND INCLUSIVE 7%(3686) MONTE 
CARLO SIMULATION 


BESIII is a general-purpose detector at the double- 
ring e*e^ collider BEPCII and is used for the study 
of physics in the 7-charm energy region [18]. It has 
a geometrical acceptance of 93% of 4» solid angle and 
consists of four main subsystems: a helium-based multi- 
layer drift chamber (MDC), a plastic scintillator time-of- 
flight (TOF) system, a CsI(TI) electromagnetic calorime- 
ter (EMC) and a resistive plate muon chamber system. 
The first three sub-detectors are enclosed in a supercon- 
ducting solenoidal magnet with a 1.0 T magnetic field. 
More details of the detector are described in Ref. [19]. 

MC simulations of the full detector are used to de- 
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termine detection efficiency and to understand potential 
backgrounds. The GEANT4-based simulation soft- 
ware, BESIII Object Oriented Simulation [1], contains 
the detector geometry and material description, the de- 
tector response and signal digitization models, as well as 
records of the detector running conditions and perfor- 
mance. Effects of initial state radiation (ISR) are taken 
into account with the MC event generator KKMC [23 [23], 
and final state radiation (FSR) effects are included in the 
simulation by using PHOTOS Bd. Particle decays are 
simulated with EVTGEN for the known decay modes 
with branching fractions set to the world average M and 
with the LUNDCHARM model for the remaining un- 
known decays. 

Angular distributions of the cascade Æ1 transitions 
(3686) — vyx«; follow the formulas in Refs. Bg, 
while the cos@ distributions for xo; — yJ/w are gen- 
erated according to phase space distributions. The Xes 
are simulated with Breit- Wigner line shapes. To account 
for the £1 transitions for (3686) > yxcj, xc > YI /v, 
MC events will be weighted as described in Section 


III. EVENT SELECTION 
A. Inclusive 1»(3686) — yX Event Selection 


We start by describing the selection procedure for 
w(3686) event candidates. To minimize systematic uncer- 
tainties from selection requirements, the (3686) event 
selection criteria, which are used for both data and the 
MC sample, are fairly loose. 

Charged tracks must be in the active region of the 
MDC with |cos0| « 0.93, where 0 is the polar angle 
of the track, and have V, < 2 cm and |V;| < 10 cm, 
where V, is the distance of the point of closest approach 
of the track to the beam line in the plane perpendicular 
to the beam line and |V;| is the distance to the point 
of closest approach from the interaction point along the 
beam direction. In addition, p « 2.0 GeV/c is required 
to eliminate misreconstructed tracks, where p is the track 
momentum. 

Photon candidates are reconstructed from clusters of 
energy in the EMC that are separated from the extrap- 
olated positions of any charged tracks by more than 
10 standard deviations and have reconstructed energy 
E, > 25 MeV in the EMC barrel (|cos6,| < 0.80) or 
> 50 MeV in the EMC end-caps (0.86 < | cos0,| < 0.92), 
where E} is the photon energy and 6., is the polar angle of 
the photon. The energy deposited in nearby TOF coun- 
ters is included in EMC measurements to improve the 
reconstruction efficiency and energy resolution. Photons 
in the region between the barrel and end-caps are poorly 
reconstructed and are not used. In addition, E, < 2.0 
GeV is required to eliminate misreconstructed photons. 
'The timing of the shower is required to be no later than 
700 ns after the reconstructed event start time to sup- 
press electronic noise and energy deposits unrelated to 


the event. 

To help in the selection of good (3686) candidates, 
events must have Na, > 0, where Nen is the number 
of charged tracks, and Ey = Een + Eneu > 0.22 Fem, 
where Ey, is the visible energy of the event, Een is the 
total energy of the charged particles assuming them to be 
pions, Eneu is the total energy of the photons in the event, 
and Ecm is the center of mass (CM) energy. To remove 
beam background related showers in the EMC and to 
demand at least one photon candidate in order to select 
inclusive (3686) — yX events, we require 0 < N, < 17, 
where N, is the number of photons. In the following, 
inclusive (3686) events and inclusive $»(3686) MC events 
will assume this selection. 


B. Non-$(3686) background 


By examining the continuum sample taken at a CM en- 
ergy of 3.65 GeV, a set of selection requirements were cho- 
sen to further remove non-7(3686) background by iden- 
tifying Bhabha events, two-photon events, ISR events, 
beam background events, electronic noise, etc. Events 
satisfying any of the following conditions will be removed: 


1. Na, < 4 and pic > 0.92Ebeam, where Epeam is the 
beam energy and the p; is the momentum of any 
charged track in the event. 


2. Na, « 4 and (EgMC)i > 0.9Epeam, where (EFemc)i 
is the deposited energy of any charged or neutral 
track in the EMC. 


3. Na, < 4 and Er < 0.15 E;4, where Feal is the 
total deposited energy (charged and neutral) in the 
EMC. 


4. Na, = 1 and (Eg, + Eneu) < 0.35 Eem. 


5. (Pan + (Pz)neu)|€ > 0.743 Epeam;, where (Pz )ch 
and (Pz)neu are the sums of the momenta of the 
charged and neutral tracks in the z direction. 


CLEO in Ref. [d] used a similar selection in their analysis. 


C. m candidate selection 


The invariant mass distribution of all yy combinations 
has a clear peak from 7° — yy decay. To reduce back- 
ground under the radiative transition peaks, photons in 
7°’s will be removed from the inclusive photon energy 
distributions. To reduce the loss of good radiative tran- 
sition photons due to accidental mis-combinations under 
the 7° peak, the requirements for a «? candidate are 
rather strict. 

Photons in 7° candidates must have 6 > 14 degrees, 
where ó is the angle between the photon and the clos- 
est charged track in the event, and the lateral shower 
profile must be consistent with that of a single photon. 
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The 7° candidates must have at least one photon in the 
EMC barrel; a one-constraint kinematic fit to the nomi- 
nal 1? mass with a x? « 200; and 0.12 « M,, « 0.145 
GeV/c?, where Mj, is the yy invariant mass. In ad- 
dition, |cos0*| < 0.84 is required for a 7° candidate, 
where 0* is the angle of a photon in the 7° rest frame 
with respect to the 7° line of flight. Real 1? mesons de- 
cay isotropically, and their decay angular distribution is 
flat. However 7° candidates that originate from a wrong 
photon combination do not have a flat distribution and 
peak near | cos 0*| = 1. 


D. Comparison of inclusive 7(3686) data and the 
MC sample 


Since efficiencies and backgrounds depend on the accu- 
racy of the MC simulation, it is important to validate the 
simulation by comparing the inclusive (3686) MC with 
on-peak data minus continuum data. In the following, 
data will refer to on-peak data minus scaled continuum 
data, where the scale factor of 3.677 accounts for the dif- 
ference in energy and luminosity between the two data 
sets a. In general, data distributions compare well with 
the inclusive MC distributions, except for those involving 
ws. To improve the agreement, each MC event is given 
a weight determined by the number of 7°s, N,o, in the 
event. For events with No corresponding to bin i of the 

data 
N,o distribution, wro = —— 

In Fig.[I]representative charged track distributions, (a) 
Nen, (b) Vz, (c) p, and (d) Egmc, are shown. Here and 
for the distributions of Figs. [2] and [8] data, unweighted 
MC, and weighted MC distributions are shown. Photon 
distributions, (a) Ny, (b) 04, (c) 6, and (d) My, of all 
yy combinations, are shown in Fig. |2} The agreement is 
acceptable for the charged distributions with or without 
weighting. For photons, the agreement for the 7° peak 
in the M,, distribution (Fig.[2] (d)) is improved with the 
weighted MC distribution, while the agreement for the 
other distributions is neither better or worse. 

Representative 7° candidate (see Section [IIT C] distri- 
butions, (a) the number of 7°s (N40), (b) the yy invariant 
mass (M.,4) made without the 7° mass selection require- 
ment, (c) | cos0*|, and (d) momentum (Po), are shown 
in Fig. B] The agreement is improved for the weighted 
sample, and in the following, the inclusive MC distribu- 
tions will be weighted by wo. 


IV. INCLUSIVE PHOTON ENERGY 
DISTRIBUTIONS 


Inclusive photon energy distributions are obtained us- 
ing the following selection requirements. First, the event 
must satisfy the inclusive 7)(3686) selection requirements, 
as described in Section [IA] and not be a non-/(3686) 
background event, as defined in Section B] a rtr J/w 
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FIG. 1. The distributions are (a) Nen, (b) Vz, (c) p, and 
(d) Ezmc. Data are represented by dots, and the MC sam- 
ple by the red and shaded histograms for the weighted and 
unweighted MC events, respectively. 


event, or a 1949 J/1» event. The +x J/w events are 
selected with the following requirements. There are 
two oppositely charged pions with momenta pr < 0.45 
GeV/c, and the mass recoiling from the ttm system, 
RM*-, must satisfy 3.00 < RM*- < 3.11 GeV/c’. 
The 292? J/4» events must have two 7°s with p, < 0.45 
GeV/c, and the mass recoiling from the nr’ system, 
RM", must satisfy 3.085 < RM? < 3.12 GeV/c?. 
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FIG. 2. The distributions are (a) Ny, (b) 64, (c) ô, and (d) 
Myy of all yy combinations. Here 0, is the polar angle of 
the photon. Data are represented by dots, and the MC sam- 
ple by the red and shaded histograms for the weighted and 
unweighted MC events, respectively. 


The photon must be in the EMC barrel. This re- 
quirement is used because the energy resolution is better 
for barrel photons, and there are fewer noise photons. 
The photon must satisfy the requirement of ô < 14 de- 
grees (see Section and not be part of a 7° can- 
didate. In Fig. [4] (a) and (b), inclusive photon energy 
distributions after the above selection requirements are 
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FIG. 3. The distributions of 7° candidates are (a) No, (b) 
My made without the 7° mass selection requirement, (c) 
| cos 0*|, (d) p,o. Data are represented by dots, and the MC 


sample by the red and shaded histograms for the weighted 
and unweighted MC events, respectively. 


shown for data and inclusive MC events, respectively. 
'The peaks from left to right in each distribution corre- 
spond to (3686) > xc, YXel YXco; Xe1 > Y/Y, 
and x2 > yJ/w. The very small peak at around 0.65 
GeV is from the (3686) — yne transition. Other small 
peaks not seen in the spectra but considered in the fit 


are J/Y > ye and Xeo > yd /W. 
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FIG. 4. Inclusive photon energy distributions for (a) data and 
(b) inclusive MC events, where the shaded region in (b) has 
the radiative photons removed. Peaks from left to right are 
w(3686) — YXc2, YXe1, YXco and xc1 and xe — yJ/w. The 


Xco — yJ/w peak is not visible. The very small peak around 
0.65 GeV is (3686) > yne. 


The inclusive 1/(3686) MC sample is used to obtain the 
signal shapes for charmonium transitions and the shape 
of the major component of the background under the 
signal peaks. The signal shape for each transition is ob- 
tained by matching the radiative photon at the gener- 
ator level with one of the photons reconstructed in the 
EMC. The requirement, which has an efficiency greater 
than 99%, is that the angle between the radiative pho- 
ton and the reconstructed photon in the EMC must be 
less than 0.08 radians. No requirement on the energy is 
used to allow obtaining the tails of the energy distribu- 
tion. The signal shapes are shown in Fig. The three 
large peaks from left to right in Fig. [5] (a) correspond to 
the (3686) — xc», YNe1, and yXco transitions. The 
very small peak around 0.65 GeV is the (3686) — Yne 
transition. The peaks in Fig. 5] (b) from left to right cor- 
respond to the y«j — yJ/w transitions for J = 0, 1, and 
2, where the y.o — yJ/w peak at around 0.3 GeV is very 
small. 

'The background component is obtained from the simu- 
lated inclusive photon energy distribution after all selec- 
tion requirements but with energy deposits from radia- 
tive photons for charmonium radiative transition events 
((3686) + vyx«j, v(3686) > We, xc; > *J/», and 
J/ > Yne) removed. Note that this distribution, shown 
as the shaded region in Fig. Ø] (b), has a complicated 
shape. This distribution will be used to describe part 
of the background under the signal peaks in fitting the 
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FIG. 5. Photon energy line shapes from inclusive MC. (a) 
Peaks from left to right are (3686) > *xc2, YXe1, and *xco. 
(b) Peaks from left to right are xco, xci, and xe2 > yJ/w. 


data and MC inclusive photon energy distributions, as 
described in Section 

The E1 transition is expected to have an energy de- 
pendence of E2, where E, is the energy of the radiative 
photon in the CM of the parent particle [29]. 'To account 
for the £1 transitions for (3686) > YXeJ, Xe3 > yJ/ b, 
a weight (trans) is calculated for each MC event using 
the radiative photon CM energy. For (3686) > «xc; 
events with no subsequent xy«; — yJ/w decay, the 


weights are given by ( = 21-)?. where Ey, for each de- 
710 


cay is the radiative photon CM energy and E410 is the 
aus E? -M2 

most probable transition energy (E10 = —5-5.22-). For 

(3686) — Y1Xes;Xer — "yaJ/ events, the weights are 


calculated according to (Fa) (2253, where E2 is the 
710 20 

energy of the daughter radiative photon in the rest frame 

of the mother particle and E~, is its most probable en- 

ergy. The overall event weight is the product of both 


weights (wro x Utrans)- 


V. (3686) + 4x.; EXCLUSIVE EVENT 
SELECTION AND PHOTON ENERGY 
DISTRIBUTIONS 


In order to constrain the final (3686) > yx-7 signal 
shapes in fitting inclusive photon energy distributions, 
clean energy spectra from (3686) > YVe7,;Xer — eX 
clusive events will be used. To fit the (3686) > YXeJ 
peaks of data, exclusive event samples are selected from 
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(3686) data. To fit the MC v(3686) > vyx.; peaks, 
exclusive samples are generated, as described below. Ex- 
clusive events must satisfy the same requirements as in- 
clusive events when constructing photon energy distribu- 
tions. 


A. (3686) — yxcz exclusive event selection 


The exclusive (3686) + *x.; photon energy distri- 
bution is the sum of (3686) > vyx«j, xc; — 2 and 4 
charged track events. 

Common requirements The number of good pho- 
tons must be greater than zero and less that 17. The 
photon with the minimum recoil, which is the angle be- 
tween the photon and the momentum recoiling from the 
two (four) charged tracks, is selected as the radiative pho- 
ton, and 0,,.,;; must satisfy 0,..,;; < 0.2 rad. Also re- 
quired are | cos, aq. .,| < 0.75, where 044. ., is the polar 
angle of the radiative photon, and 3.3 < M»(4, < 3.62 
GeV /c?, where M»(4)« is the invariant mass of the two 
(four) charged tracks. 

Specific requirements for (3686) — 
YXeJ>XeJ — 2 charged tracks We require one 
positively and one negatively charged track. Particle 
identification probabilities are determined using dE/dx 
information from the MDC and time of flight informa- 
tion from the TOF system, and both tracks are required 
to be either kaons (Prob(K) > Prob(m)) or pions 
(Prob(z) > Prob(K)). We also require |cos0| < 0.85 
for both charged tracks, where 0 is the polar angle, the 
momentum of each track be less than 1.4 GeV/c, and 
the momentum of one track be larger than 0.5 GeV/c. 

Specific requirements for (3686) — 
C"XcJ:XcJ — 4 charged tracks We require two 
positive and two negative tracks and |Xp;| < 0.04 
GeV/c, where |p-| is the sum of the momenta of the 
charged tracks and neutral clusters in the z direction. 
ISR events tend to have large |Xp.|. Also the mass 
recoiling from the two low momentum tracks is required 
to be less than 3.05 GeV/c? to veto (3686) > amJ/w 
background. 


B. v(3686) > yxc; exclusive MC sample 


Here, exclusive yez — two and four pion and kaon 
events are generated with EVTGEN [23], and the gener- 
ated events are selected using the selection criteria de- 
scribed in Section Events are weighted by wirans 
using the generated energy of the radiative photon. 


VI. FITTING THE INCLUSIVE PHOTON 
ENERGY DISTRIBUTION 


The numbers of (3686) — yx.; events and x. > 
*J/' events are obtained by fitting the inclusive photon 


energy distributions for data. The efficiencies are ob- 
tained from the fit results for the inclusive w(3686) MC 
events. 


To fit the (3686) + vyx.; signal peaks of data, the 
MC signal shapes, described in Section are convolved 
with asymmetric Gaussians to account for the difference 
in resolution between MC and data, where the param- 
eters of the Gaussians are determined by the fit. The 
broad x«i and xc — ^4J/w peaks are described well 
by just the MC shapes. Also included in the fit are 
Xco > *J/? and J/ > yne. The background distribu- 
tion is the inclusive MC photon energy distribution with 
energy deposits from radiative photons removed com- 
bined with a second order Chebychev polynomial func- 
tion. 


To constrain further the (3686) — YXez signal 
shapes, a simultaneous fit to inclusive (see Section [IV] 
and exclusive photon energy distributions (see Sec- 
tion [VA] is done in the energy range from 0.08 to 0.35 
GeV. The parameters of the asymmetric Gaussians are 
the same for the inclusive and exclusive fits. However, all 
signal shapes are allowed to shift independently in energy 
for the two distributions. The exclusive background dis- 
tribution is determined in a similar way as the inclusive 
photon background distribution but using the exclusive 
event selection on the (3686) MC event sample. 


Shown in Fig. [6]is the simultaneous fit of data for the 
region 0.08 < E, < 0.5 GeV for the inclusive photon en- 
ergy distribution and the region 0.08 < E, < 0.35 GeV 
for the exclusive photon energy distribution. The fit to 
the inclusive photon energy distribution and the corre- 
sponding pull distribution are shown in the top set of 
plots. The bottom set of plots are those for the exclusive 
photon energy distribution. The pull distributions are 
reasonable, except in the vicinity of the (3686) > yxe1 
and vyxe» peaks. The chi-squares per degree of freedom 
(ndf) axe 3.5 and 2.7 for the inclusive and exclusive dis- 
tribution fits, respectively. The chi-square is determined 
using X? = X;((n; — n)/o;)?, where ni, nf, and o; are 
the number of data events in bin i, the result of the fit at 
bin z, and the statistical uncertainty of n;, respectively, 
and the sum is over all histogram bins. 


A fit is also done to the MC inclusive energy distribu- 
tion. The MC shapes are used without convolved asym- 
metric Gaussians for the (3686) > YXez peaks. Since 
only MC shapes are used, it is not useful to do a si- 
multaneous fit as there are no common parameters to 
be determined in such a fit. The fit matches the inclu- 
sive photon energy distribution almost perfectly with a 
chisquare close to zero. This is not unexpected since the 
signal and background shapes come from the MC and 
when combined reconstruct the MC distribution. 
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TABLE II. Branching fraction results. The indicated uncertainties are statistical only. 


Branching Fraction 


B(1(3686) > yxco) 

B(4(3686) — yxe1) 

B(V(3686) > 7x2) 

B(v(3686) — xc0) x B(xco — 4J/) 
B(4(3686)  yxaà) x B(xa > yJ/v) 
B(w(3686) > yxe2) x B(xea + yJ/) 
B(xco + yJ/V) 


= 
eo 


TTPO [ETET [EAS 


Events/0.5 MeV 


Events/1 MeV 


FIG. 6. Simultaneous fits to the photon energy distribu- 
tions of data. (Top set) Inclusive distribution fit and corre- 
sponding pulls, and (Bottom set) exclusive distribution fit and 
pull distribution. Peaks from left to right in the top set are 
10(3686) > YxXe2, YXe1, and yXco and xc and xe» > yJ/wv. 
The xco — *J/w peak is not visible. The smooth curves 
in the two plots are the fit results. The dashed-dotted and 
dashed curves in the top plot are the background distribu- 
tion from the inclusive 1/(3686) MC with radiative photons 
removed and the total background, respectively. The back- 
ground in the exclusive fit plot is not visible. 


0.0123 + 0.0081 
1.8881 + 0.0053 3.442 + 0.010 
0.9828 + 0.0041 


Events (x109) |Efficiency| Branching 
Fraction (96) 


- 0.0068 
- 0.0054 
- 0.0055 


9.389 + 0.014 
9.905 + 0.011 
9.621 + 0.013 


0.024 + 0.015 


1.793 + 0.008 


0.25 + 0.16 
34.75 + 0.11 
18.64 + 0.08 


VII. BRANCHING FRACTION 
DETERMINATIONS 


The branching fractions are calculated using the fol- 
lowing equations 


N, 
B(Y(3686) > yxer) = MM. 
€15(3686)—5x«, X 1 'w(3686) 


where B((3686) — -x.;) is the branching fraction 
of (3686) — Yxes, Nw(3686)54x., i$ the number of 
events in data from the fit, €4(3686).,45,, is the effi- 
ciency determined from MC, and (3686) is the number 
of (3686) events [13. 'The product branching fraction 
for (3686) > YXcJ, Xe — *J/v is given by 


B(v(3686) > vyxcj) x B(xcs > ?J/w) 


= Ny Jf) (2) 


€y,524J/o X Ny(3686) 


where N,.,. ,,7/, is the number of xo; — ^J/» events 
in data and e,,, ,,7/, is the efficiency. From Eq. (1) 
and Eq. (2), we obtain the branching fraction for y.7 > 
yJ/w, which is given by 


B(xer > 4I/4) 

B(Y(3686) > Xes) x B(xes > S/W) 
B(4(3686) > Xes) 

€(3686) 7x07 X NyesyI/v 


(3) 


Exes oy I/b X Nu (3686) 7x04 


Results are listed in Table where the uncertainties 
are statistical only. For B(Xez — yJ/W), an alternative 
parametrization in terms of N,(3686).,,4,, and the ra- 
tio Ny, +47/v/Nuw(3686)>7x., has been tried because of 
the possible correlation between the numerator and de- 
nominator of Eq. (3), but the difference with the original 
result is small and will be neglected since it is much less 
than the systematic uncertainties that will be discussed 
below. 
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VIII. SYSTEMATIC UNCERTAINTIES 


Systematic uncertainties, which arise from selection re- 
quirements, fitting, photon efficiency, the uncertainty in 
the number of (3686) events, etc. are summarized in 
Table [III] For (3686) > ^x;;, they are under 4% and 
smaller than those of CLEO |9], with the largest contribu- 
tion coming from fitting the photon energy distribution. 
Details of how they are estimated are given below. 


A. Systematic uncertainties from initial (3686) 
event selection 


Initial (3686) event selection requirements are Nen > 
0, Ny < 17, and Ey > 0.22Ecm. To determine the 
systematic uncertainties associated with the Nen > 0 re- 
quirement, events without charged tracks are also ana- 
lyzed. The photon time requirement is removed for these 
events since without charged tracks, the event start time 
can not be well determined. The selection requirements 
are also changed because these events have much more 
background. Events must have total energy greater than 
1.7 GeV and at least one good neutral pion. Even so, 
there is a background from low energy photons, and even 
after subtracting continuum, the photon energy distri- 
bution for data has a large background under the sig- 
nal peaks, making fits difficult with the number of fitted 
events having large uncertainties. 

The photon energy distributions for data and MC are 
fitted. The numbers of fitted events for data and MC are 
then added with the number of fitted events with charged 
tracks, and the branching fractions are recalculated. The 
differences with the branching fractions determined with 
charged track events only are then determined and taken 
as the systematic uncertainties associated with the Nen > 
0 requirement. 

As described in Section inclusive (3686) MC 
events are weighted according to the the N,o distribution 
to give better agreement with data. According to the 
MC, the efficiency of the N, < 17 requirement, defined 
as the number of events with N, < 17 divided by the 
number of events with N, > 0, is 99.99% with weighting 
and 99.99% without weighting, while the efficiency for 
data is 99.98%. The agreement is excellent, the efficiency 
is very high, and the systematic uncertainty is negligible 
for this requirement. 

The agreement between the Ey, distribution of data 
and the inclusive 1/(3686) MC distribution is very good. 
According to the inclusive MC, the efficiency of the 
Eyis > 0.22Ecm requirement after the Nen and Ny re- 
quirements is 99.76%. The mean and root-mean-squared 
values of the MC (data) are 3.004 (2.991) GeV and 0.561 
(0.579) GeV, respectively. If the MC distribution is 
shifted down by 13 MeV relative to the data, the loss 
of events due to the Ey;, requirement corresponds to an 
inefficiency of 0.17%, and this will be taken as the sys- 
tematic uncertainty for the Eyis requirement. 


10 


B. Systematic uncertainties from inclusive photon 
selection 


Further selection criteria are used before including 
photons into the photon energy distributions which are 
used for fitting. Photon selection requirements include 
ô > 14°, removal of non-7)(3686) background events, re- 
moval of m7zJ/y events, and removal of photons which 
can be part of a 7°. 

ô > 14? and v(3686) background removal sys- 
tematic uncertainties: 

To determine the systematic uncertainties for the first 
two requirements, they are removed from the selection 
process, and the branching fraction results obtained are 
compared to those with the requirements. Removing the 
6 requirement changes the inclusive photon energy back- 
ground distribution of the inclusive MC, as well as the 
inclusive photon energy distribution of data. The differ- 
ences of the branching fraction results are taken as the 
systematic uncertainties for each of the requirements. 

atn—J/w event removal systematic uncer- 
tainty: 

The distribution of mass recoiling from the v * «^ sys- 
tem, RM'-, for events passing the non-v(3686) veto 
and the (3686) + ntn J/w selection, but without the 
recoil mass requirement in Section has a clear J/y 
peak from 4(3686) — «*« J/w. Events with RM- 
satisfying 3.09 < RM*- < 3.11 GeV/c? will be removed 
from further consideration. However, there are «^-^ 
mis-combinations underneath the peak in the J/ region. 
'To estimate the probability that a good radiative photon 
event may be vetoed accidentally (or the efficiency with 
which it will pass this veto requirement), the sideband 
regions, defined as 3.07 < RM*~ < 3.085 GeV/c? and 
3.115 < RM*~ < 3.13 GeV/c?, are used to estimate the 
number of mis-combinations in the signal region. Us- 
ing this veto probability, the efficiency for inclusive MC 
events to pass the (3686) > ata J/wv veto require- 
ment is found to be 93.06%. The efficiency for data 
is 92.83%, and the difference between data and inclu- 
sive MC is 0.23/93.06 — 0.2596, which we take as the 
systematic uncertainty due to the «^-^ J/ veto for all 
radiative photon processes. 

TOTO J /1b event removal systematic uncertainty: 

The approach to determine the systematic uncertainty 
for the 7°n°J/y event removal is similar to that de- 
scribed in the previous section. Using the veto probabil- 
ity obtained using sidebands, the efficiency for inclusive 
MC events to pass the /(3686) — 1929/4» veto require- 
ment is found to be 95.34%. The efficiency for data is 
95.37%, and the difference between data and inclusive 
MC is 0.03/95.35 = 0.03%, which we will take as the 
systematic uncertainty due to the n?n? J/w veto for all 
radiative photon processes. 

Systematic uncertainty for the removal of pho- 
tons which can be part of a 7°: 

As described in Section [LH C] photons that are part of 
a 1? are excluded from the inclusive photon energy distri- 
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TABLE III. Relative systematic uncertainties (%). By z is notation for B((3686) > yxe7), Bes is for B((3686) > yxe.r) x 
B(xer > yJ/w), and By; is for B(xc; — yJ/wW). Some uncertainties cancel in the determination of B44 and By2 and are left 
blank in the table. Since the fit uncertainty is so large for (3686) > yxco, xco — yJ/w, the other systematic uncertainties for 


Bpo and Byo are omitted. 


Buyo Bui 


Na > 0 

Ny <17 

Evis = 0.22 Ec 

6 > 14? 

~)(3686) background veto 
ata Jf veto 

n^n? Ji) veto 

^y not in x? 

Fitting 

MC signal shape 
Multipole correction 

| cos 0| « 0.8 

v? weight 

Continuum energy difference| 0.75 
^y efficiency 
IN.) (a686) 


Total 


bution. To estimate the systematic uncertainty for this 
requirement, the efficiency of this criterion is determined 
for data and MC events for each transition by fitting the 
photon inclusive energy distribution with and without 
the 7° removal using non-simultaneous fitting. The sys- 
tematic uncertainties are determined by the differences 
between the efficiencies for data and MC events. 


C. Fitting systematic uncertainty 


The systematic uncertainty associated with the fit pro- 
cedure is determined by comparing various fitting meth- 
ods. The fit is done with an alternative strategy, fitting 
with a non-simultaneous fit, changing the order of the 
polynomial function used from second order to first or- 
der, changing the fitting range, and fixing the number of 
events for the J/w — ^n. and (3686) > *xco, Xco > 
yJ/w to the numbers expected, and the result for each 
case is compared with our standard fit to determine the 
systematic uncertainty for that case. 

For the alternative strategy, the (3686) > yxe1 and 
^Xc2 peaks are described by asymmetric Gaussians with 
Crystal Ball tails on both sides. The other signal peaks 
and backgrounds are the same. A simultaneous fit is 
done to better constrain the asymmetric Gaussian and 
Crystal Ball tail parameters, which are common between 
the inclusive and exclusive distributions. 

For the 1/(3686) + »yx«; systematic uncertainties, the 
fitting range is changed from 0.08 — 0.5 GeV to 0.08 — 
0.35 GeV, which removes the x.; — yJ/w peaks and 
changes the number of parameters used in the fit. For 
the xe; — yJ/w systematic uncertainties, the range is 


By2 | Bro Bp Bp2 | Byo By Bye 


0.06 0.74 0.21 1.5 


0.17 0.17 
0.02 1.56 0.12 1.42 
0.51 0.11 1.25 0.26 
0.25 0.25 
0.03 0.03 
1.24 2.3 1.35 2.3 
3.10 7.22 | 861 4.43 7.27 
0.07 0.53 0.24 1.05 
0.35 3.82 0.70 3.87 
0.35 1.46 0.47 1.52 
1.09 1.73 0.47 0.13 
0.35 0.60 0.39 1.02 
10 10 
0.81 0.81 


3.54 3.57 2.83 | 869 3.84 9.09 | 861 4.92 9.05 


changed from 0.08 — 0.5 GeV to 0.2 — 0.54 GeV, which 
removes the 1/(3686) + yXcı and (3686) > yx-2 peaks 
and produces a rather large systematic uncertainty due to 
the background in the fit of data preferring a pure polyno- 
mial background in the latter case. The total systematic 
uncertainties from fitting for each branching fraction are 
determined by adding the systematic uncertainties from 
each source in quadrature. 

The signal for (3686) > YXco, Xco — YyJ/1b is very 
small and sits on the tail of (3686) > yx-o. It is there- 
fore difficult to fit this peak as indicated by the very large 
fitting systematic uncertainty for this process. 


D. MC Signal Shape 


The signal shapes used in fitting the photon energy 
distribution are determined by matching MC radiative 
photons with reconstructed photons in the EMC, where 
the angle between the photons is required to be less than 
A0 — 0.08 radians. This selection could bias the signal 
shapes used in the fitting. The systematic uncertainty as- 
sociated with this selection is determined by changing the 
A9 selection requirement to 0.04 radians. The differences 
for each decay are taken as the systematic uncertainties 
in the signal shape. 


E. Higher order multipoles for (3686) > yxc1 
and Xc2 


Angular distributions for (3686) > »x.; are gener- 
ated according to those expected for El radiative tran- 


201711.02276v1 


chinaXiv 


sitions. This approach is accurate enough for (3686) > 
^Xco; but higher order multipole contributions must be 
considered for (3686) > yxXe1 and (3686) > »yx«» de- 
cays. Also the angular distributions for xo; — yJ/w 
MC events do not agree with data. BESIII has measured 
the angular distributions for (3686) > yxer,XNer > 
yJ fib lig. and these distributions have been fitted to 
1+a cos? 0, where 0 is the laboratory polar angle, and the 
values of a have been determined. Using these values of 
a, it is possible to calculate the differences in the geomet- 
ric acceptance between data and the inclusive (3686) 
MC. The acceptance efficiency for a given value of a is 
given by the integral of 1 + a cos? 0 from cos 0 = —0.8 to 
cos = 0.8 divided by the integral between —1 and +1. 
Using the values of a that were used to generate the MC 
events and those obtained based on Ref. , the changes 
in the efficiencies are 0.61% for (3686) > yxe1 and 
0.60% for (3686) — yxe2. For xez > yJ/w (J = 1,2), 
the changes are 0.35% and 3.82%, respectively. The 
changes to the branching fractions from the changes in 
efficiencies are taken as the systematic uncertainties due 
to the higher order multipole corrections. 


F. |cos0| < 0.8 


The systematic uncertainty associated with the 
[cos0| < 0.8 requirement is determined by using the 
requirement |cos0| < 0.75 instead and by comparing 
the results with the standard requirement. This tests 
whether there are edge effects with the EMC that are 
not fully modeled by the MC simulations. 


G. Event weighting 


As described in Section [ID] MC events are weighted 
to give better agreement for the 7° distributions between 
data and MC simulation, as well as to include the E1 
transition E3 weight. The systematic uncertainty asso- 
ciated with the w,o weight is determined by turning off 
its weighting and taking the difference in results as the 
systematic uncertainties. 


H. Continuum energy difference 


Data distributions, including the inclusive photon en- 
ergy distribution for data, are defined as data minus 
scaled continuum data. While this takes into considera- 
tion the effect on the normalization of the continuum due 
to the difference in luminosity and energy, it does not 
consider the difference in the energy scale of the pho- 
tons. To determine the systematic uncertainty due to 
this effect, the photon energies of the continuum data 
were scaled by the ratio of the CM energies, 3.686/3.65, 
and the scaled distribution was subtracted from data, 
and the fitting redone. The differences with respect to 
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the standard analysis are taken as the systematic uncer- 
tainties of this effect. 


I. Other systematic uncertainties 


'The photon detection efficiency is studied utilizing the 
control samples 1/(3686) 3 rtr J/q», J/i — pn? and 
(3686) — «919 J/V with J/4 > I*I- and ps9. The 
corresponding systematic uncertainty is estimated by the 
difference of detection efficiency between data and MC 
samples, and 1% is assigned for each photon [ad]. 

The trigger efficiency is assumed to be very close 
to 10096 with negligible uncertainty, since the average 
charged particle and photon multiplicities are high. The 
number of 1/(3686) events is (106.414-0.86) x 10°, which is 
obtained by studying inclusive (3686) decays m The 
uncertainties from all above sources and the total sys- 
tematic uncertainty, obtained by adding all uncertainties 
quadratically, are listed in Table [NI] Since the fitting un- 
certainty for (3686) > YXco;, Xeo — ^yJ/w is so large, 
indicating that this fit is not very meaningful, only this 
uncertainty is listed in the table. 


IX. RESULTS 


Our results are listed in Table We also calculate 
ratios of branching fractions, where common systematic 
uncertainties cancel 


B(v (3686) > YXco)/B(Y (3686) — 7x1) 
= 0.948 + 0.002 + 0.044 


B(v (3686) > YXco)/B(Y (3686) => x2) 
= 0.976 + 0.002 + 0.040 
B(v(3686) > Xc2)/B((3686) => x1) 
= 0.971 + 0.002 + 0.040 


For comparison with some theoretical calculations, we 
also determine partial widths using our branching frac- 
tions and the world average full widths i. Table 
contains our partial width results, as well as theoretical 
predictions, reproduced from Table VI in Ref. gj. 'The 
theoretical predictions include the linear potential (LP) 
and screened potential (SP) models [s], as well as earlier 
predictions from a relativistic quark model (RQM) [33], 
non-relativistic potential and Godfrey-Isgur relativized 
potential models (NR/GI) [34], and color screened mod- 
els, calculated with zeroth order wave functions (SNRo) 
and first order relativistically corrected wave functions 
(SNR1) [B5]. 


X. SUMMARY 


Our results, CLEO measurements id, [31] [33], previ- 
ous BESIII measurements [1d], and PDG results 
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TABLE IV. Our branching fraction results, other results, and PDG compilation results. 


Branching Fraction This analysis (96) 


B(w(3686) + yxca) x B(xa — yJ/) 3.442 + 0.010 4 


B(1(3686) > yxe2) x B(xea — yJ/%) | 1.793 + 0.008 4 


Xco > *J /i) 
B(xa > yJ/A) 


Other 


9.22 + 0.11 + 0.46 [9] 
+0.11 + 0.54 [9] 
+ 0.14 x 0.61 [9] 


9.07 
9.33 


v(3686) — *xco 9.389 + 0.014 4 
v(3686) + *xa 9.905 + 0.011 4 
v(3686) + 7X2 9.621 + 0.013 4 
(3686) — *yxco) x B(x«o — yJ/w)|0.024 + 0.015 + 0.205] 0.125 4 


0.151 4 


- 0.007 + 0.013 
- 0.003 Æ 0.010 


0.158 d 


E 0.003 Æ 0.006 


3.56 + 0.03 + 0.12 [3 


3.377 d 


3.518 + 0.01 d 
1.95 + 0.02 d 


1.874 d 


E 0.009 + 0.183 


- 0.007 + 0.102 


1.996 4 


- 0.008 Æ 0.070 


- 0.120 [16] 
- 0.07 [31] 


+ 0.035 0.127 + 0.006 


3.24 + 0.07 


1.75 + 0.04 


33.9 + 1.2 


Xc2 — yJ/b) 


TABLE V. Partial widths (keV) of radiative transitions for (3686) > yJ/w and xc; — yJ/wv. Shown are our experimental 
results and predictions from a relativistic quark model (RQM) ; non-relativistic potential and Godfrey-Isgur relativized 
potential models (NR/GI) [34]; color screened models [33], calculated with zeroth order wave functions (SNRo) and first order 
relativistically corrected wave functions (SNR); and linear potential (LP) and screened potential models (SP) ER The zi 
predictions include only F1 transition calculations, while the za results include higher order multipole corrections. 


Initial | Final Imi (keV) Tem (keV) 
state |state |RQM [33] NR/GI [34] SNRo,, [35] LP [8] SP [8] | LP [8] SP [8] This analysis 
v(3686) 63/26 74/25 27 22 32 269-18 
54/29 62/36 45 42 45 283-19 
38/24 43/34 36 38 46 275-17 
Xc0 152/114 167/117 141 172 179 
Xe 314/239 354/244 269 306 319 306 + 23 
Xc2 424/313 473/309 327 284 292 363 + 41 


are listed in Table [V] Our (3686) — yyey branching 
fractions are the most precise. The branching fractions 
for (3686) > yxez7 agree with CLEO within one stan- 
dard deviation, except for (3686) — yxXe1 which dif- 
fers by 1.3 standard deviations. The product branch- 
ing fractions B(v(3686) > 7x1) x B(xa — yJ/1) and 
B(w(3686) > yxXe2) x B(xea — yJ/Y) agree with the 
previous BESIII measurements. Because of the difficulty 
in fitting (3686) > yxe0,Xe0 — J/w, our product 
branching fraction has a very large systematic error com- 
pared with those using exclusive decays. 

Partial widths are shown in Table[V] For comparison 
with models, experimental results have become accurate 
enough (partly due to this measurement) to become sen- 
sitive to fine details of the potentials, e.g. relativistic 
effects, screening effects, and higher partial waves. 
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